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ABSTRACT 

Optimal  control  of  the  deformable  mirror  in  a  High  Energy  Laser  was 
used  to  reduce  the  effects  of  atmospheric  turbulence.  The  performance 
characteristics  of  an  optical  detector  with  one  or  four  sub-apertures  were 
estimated.  These  estimated  performance  characteristics  were  then  used  to 
determine  the  regime  where  four  sub-apertures  gave  better  performance 
than  one  sub-aperture  for  no  time  delay  and  for  a  one  millisecond  time 
delay. 
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I.  INTRODUCTION: 


Optimal  control  of  the  deformable  mirror  in  a  High  Energy 
Laser  (HEL)  (see  Figure  1)  is  used  to  offset  the  turbulence  in  the 
atmosphere.  This  research  effort  seeks  the  best  control  system  for  the 
HEL  as  a  function  of  the  measurement  statistics  and  the  number  of  detector 
sub-apertures.  The  number  of  sub-apertures  is  important  because  more 
sub-apertures  allow  the  estimation  of  higher  order  optical  effects  and 
the  potential  for  better  control.  However,  increasing  the  number  of 
sub-aperture  increases  the  measurement  variance. 

The  primary  tool  used  in  this  determination  was  computer  modeling. 
Kleinman’s  routines  (References  1  through  5)  were  used  to  determine  the 
optimal  continuous  control  and  the  steady-state  gain  of  a  continuous 
Kalman  estimator.  Kleinman*s  routines  were  also  used  to  find  the  overall 
system  covariance  matrix. 

II.  OBJECTIVES  OF  THE  RESEARCH  EFFORT: 

The  objective  of  this  project  was  to  determine  the  system 
for  the  deformable  mirror  which  minimized  the  effects  of  atmospheric 
turbulence.  The  reduction  in  atmospheric  turbulence  effects  was  found 
as  a  function  of  measurement  statistics  and  number  of  sub-apertures. 

The  reduction  was  found  for  the  first  five  optical  modes  with  time  delays 
of  0  and  1  msec.  A  basic  assumption  was  that  each  optical  mode  could  be 
estimated  separately  and  linearly. 

III.  ATMOSPHERIC  TURBULENCE: 

The  distortion  induced  in  the  laser  beam  by  atmospheric  turbulence 
can  be  described  in  terms  of  its  effect  on  the  various  optical  modes  (tilt, 
focus,  astigmatism,  etc.).  The  effect  on  each  optical  mode  can  be  modeled 
by  a  three  state  variable  model  driven  by  a  Gaussian  random  variable 
(References  6  and  7). 


Simplified  HEL  Control  System 


Let  *  random  effect  on  optical  mode 

X2  =  turbulence  model  state 

=  turbulence  model  state 
Then  the  state  model  is 


-  a  Xj  +  a  X2 

(1) 

^2 

-  b  X2  +  b  x^ 

(2) 

II 

-  C  X3  +  v^ 

(3) 

where  v^  is  white  noise  with  covariance 

E  {v^  (t)  v^.  (t  +  T)}  =  k6  (T) 

(4) 

Values  of  k,  a, 

b,  and  c  for  the  first  five  optical 

modes  are 

given  in 

Table  I  (Reference  7). 

Table  I.  Parameters  for 

Turbulence 

Model 

Mode 

Name  k 

a 

b 

c 

X 

Tilt  2.6  X  10“^ 

47.10 

2200 

2200 

y 

Tip  2.4  X  10  ^ 

47.10 

2200 

2200 

2.2 

X  +  y  -  R  /2 

Focus  3.3  X  10  ^ 

94.20 

3350 

3350 

2  2 

X  -  y 

Astig-1  5,4  X  10  ® 

41.80 

2576 

2576 

xy  (exact) 

Astig-2  6.0  X  10 

4.710 

73.30 

2932 

xy  (approx.)  5.65  x  10  ® 

73.30 

2932 

3  X  10** 

NOTE:  Mode  xy  (exact)  has  a  zero  at  .0523  radians/sec 

The  steady-state  output  covariance  of  the  system  described  by  equations 
(1)  through  (4)  may  be  found  by  solving  the  linear  variance  equation 
(reference  9) 

A  X+  a''’x  +  Q  =  0  (5) 

where 

A  is  the  3  by  3  system  matrix  from  equations  1  through  3 

Q  is  a  3  by  3  matrix  with  q^o  =  k 
the  only  non-zero  term 

T 

X  is  the  covariance  matrix  of  [x^  X2  x^] 


The  resultant  entries  of  x  are 


V  -  X  =  X  -  ab(a->-b+c)k  (6) 

n  ■  12  *21  "  2(a+b)(a+c)(b+c)c 

X  =  X-  =  ab  k  (7) 

2(bc) (a+c)c 


*22  ~  *23  ■  *32  ”  2(b+c)  c 


(8) 


’  (9) 

JJ  -  2c 

To  avoid  a  special  program  for  mode  xy  (astigmatism-2),  mode  xy  was 

modelled  by  the  approximate  model  shown  as  the  last  line  of  Table  I.  This 

approximate  model  resulted  from  two  approximations.  First,  the  zero 

at  ,0523  was  cancelled  by  the  pole  at  4.71.  (This  cancellation  multiplied 

k  by  4. 71/. 0523).  Secondly,  the  linear  variance  equation  for  the  resulting 

2  pole  system  was  solved  to  yield 

-  a  (10) 

11  2(a+b)b 


Finally,  an  effective  k  was  chosen  so  that  the  variance  of  from  (6)  and 
the  variance  from  (10)  were  identical  when  a  remote  third  pole  (c=10b) 
was  included. 

The  effect  of  atmospheric  turbulence  upon  the  intensity  at  the  target 
is  a  function  of  Zernike  residual  errors  (reference  11  and  13).  These 
residual  errors  and  differences  are  tabulated  in  Table  II. 


Table  II.  Zernike  Residual  Errors  and  Strengths 


Mode 

A. (=o  2) 

Difference 

1.03 

(Strength  of  aberration) 

Tilt 

0.582 

.448 

Tip 

0.13A 

.448 

Focus 

0.111 

.023 

Astig-1 

0.088 

.023 

Astig-2 

.023 

5/3 

All  quantities  must  be  multiplied  by  (D/r^) 

4 


The  normalized,  on-axis  for  field  irradiance,  is  usually  found 

from  the  Strehl  approximation 

^rel  = 

The  Strehl  approximation  is  not  valid  for  large  values  of  a^.  Thus,  the 
figure  of  merit  used  in  this  research  effort  is  the  normalized  phase 
variance  where  =  1.03  (D/r^)^'  .  For  example,  if  tilt 

and  tip  could  be  completely  corrected  (or  removed)  the  normalized 
variance  would  be  0.134/1.03. 


IV.  MIRROR  AND  CONTROL  SYSTEM  MODEL: 

The  mirror  tilt  and  control  system  are  represented  by  three  state 

variables . 


Let  x^ 

=  mirror  tilt 

(11) 

^5 

=  mirror  tilt  rate 

(12) 

^6 

Finally, 

=  integral  of  error  (x^  -  x^) 

the  mirror  drive  is  a  force  driver  such  that 

(13) 

^4 

=  X5  =  u 

(14) 

The  complete  tilt  system  model  is 


j_  -a  a  0  0 

0  -b  b  0 

0  0  -c  0 

0  0  0  0 

0  0  0  0 

-  1  0  0-1 


0 

0 

0 

1 

0 

0 


0^ 

r^ii 

0- 

-0^ 

0 

X2 

0 

0 

0 

Xo 

+ 

0 

u  + 

v^ 

3 

t 

0 

^4 

0 

0 

0 

^5 

1 

0 

0- 

Lx,  J 

_o. 

.0  - 

(15) 


The  covariance  matrix  for  the  driving  noise  is  all  zeros  except  for 

*^33  ~  ^  (16) 
The  remaining  four  optical  modes  will  be  modeled  by  the  same  general 

equation.  The  different  values  from  Table  I  will  be  used  in  equations  (15) 
and  (16). 


V.  SYSTEM  OBSERVATION  MATRIX: 

The  observation  system  will  be  represented  hy  y  =  M  x  +  v 
The  quantity  yj  is  measured  by  the  detector. 

Yl  =  *1  -  +  Vi  (17) 

Observability  requires  that  also  be  measured 

73  =  ^6  +  V3  (18) 

The  programs  developed  for  the  research  (reference  14)  also  permitted 
the  measurement  of 

V2  =  *4  +  ^2 

However,  a  knowledge  of  x^  did  not  improve  the  system  performance 
and  the  coefficients  associated  with  this  measurement  were  later  set 
to  zero. 
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VI.  OPTIMAL  CONTROL  SYSTEM  DESIGN: 

The  classical  infinite  time  optimal  control  problem  minimizes  a 
cost  function  (Reference  8) . 

V  =  QjX  +  RjU)  dt  (20) 

The  minimization  of  V  produces  the  control  in  terms  of  ,  the  solution 
to  a  Ricatti  equation 

PjA  +  A^Pj  -  Pj  b  Pj  +  Qj  =  0  (21) 

(Kleinman's  subroutine  MRIC  solves  this  equation.) 

The  optimal  control  is 

u*  =  -  Rj‘^  Pj  (22) 

The  optimal  control  can  be  required  to  have  a  prescribed  stability  by 
replacing  A  by  (odJ  +  A). 


VII.  KALMAN  STATE  ESTIMATOR: 

If  the  basic  plant  is  described  by  (reference  8) 

x  =  Ax  +  Bu  +  v 
-  -  ”  -x 

y  =  M  X  +  V 
i  -  -m 

The  Kalman  estimator  is 

X  =  Ax  +  Bu  +  k  (Mx  -y) 

Desired  control  law 

u  =  L  x 

Actual  control  law 


u  =  L  X 
-  -  e 


Defining 


e  =  X  X 

~  —  —  e 

yields  a  new  system  with  doubled  dimension 


-  X  - 

rA  +  B  L 

-  B  L-i 

-X  - 

r  V  1 

• 

e 

= 

0 

A  +  k  M 
e  J 

e 

+ 

-X 

V  +  k  V 

L  -X  e  -m-i 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 
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(30) 


The  matrix  is  the  Kalman  gain.  It  is  given  by 

k  =  -  P.  M  ^  R 

e  2  m 


where  h  is  the  covariance  of  v  . 
m  -m 


(the  covariance  of  e)  is  the  solution  to  (Reference  8) 
P..  A  +  A  Po  -  P.,  M  R  M  P..  +  I 


=  0 


(31) 


Q  is  the  covariance  of  v  . 


R  is  the  covariance  of  v  . 
m  -m 


The  input  covariance  for  equation  (29)  is 


cov 


1 

=  E 

> 

u 

> 

> 

> 

'X  -X  -X  ~ffl  e  1 

V  +  k  V 

V  +  k  V  * 

L-x  e  -ffl  J 

L-x  e  Hm 

Thus , 


cov 


rv 

1  V  + 

k  V 

[Qn 

Q 

Q  + 

tl 

k  R  k  ^  J 

L-x 

e  -m-* 

L^n 

^n 

e  m  e 

Qi2 


(32) 


The  (2,2)  term  can  be  further  simplified  by  substituting  equation  (30). 

2^ 

4^  D  13 

2 


k.  R_  k_"  =  PoM"  R_  "  R_  (R„^^)*^  M  P2 


e  m  ,  m  '  m 

=  P2(m’^  R^’'  M)  P, 


Thus,  the  (2,2)  term  becomes 


Q  +  k  R  k  =  Q  +  P„  (M*^  R  M)  P„ 
^n  erne  ^n  2  m  ^2 


(33) 


The  steady-state  output  covariance  for  the  double  dimension  system 
described  by  equation  (29)  is  the  solution  of  the  linear  variance  equation 

(34 


^12  ^12  ^12  ^12  ^12  ®  ^  ^ 
where  ?^2  covariance  of  [x^  e^] 


A^2  is  the  double  dimension  system  matrix 
Qj2  is  the  input  covariance  from  (32) 

Kleinman  (reference  10)  modified  the  system  described  above  by 
including  the  effects  of  time  delay  in  the  observation  system. 


y  (t)  =  M  X  (t  -  T)  +  v  (t  -  T) 


(35 


Given  this  observation  model  Kleinman  obtains  the  following  expression 
for  the  covariance  of  x. 


^  AT^  A^T. 
p  =*  e  ^ 

X  ^ 


L 


*r« 

/o 


Me^  ^  dt 

2  m 


(36) 


where  A  =  A  +  BL  (term  from  (29)) 

(Kleinman' s  routine  INTEG  performs  numerical  integrations  of  this  type). 


F 


VIII  DETECTOR  MODELS: 

Martin  (reference  11)  presents  this  result  for  the  variance  of  an  upgraded 
shearing  interferometer  measurement 


o  9 

TT" 

2 

^"det  1 

'  VfovV  I 
( ;  1 

SI 

sin"^  j 

rTT  a 

o 

‘^FOV 

CiDf 

1 

(37) 


where  =  4. 17  x  10  ^ 

=  number  of  sub-apertures 

=  4  X  10  ^  radians 

“fOV  =  2  X  10“^  M 

Of  is  shown  in  Table  III. 
o 

Substituting  values  into  equation  (37)  produces 

a  2  =  f  (38) 

sin^Ig/M] 

With  f  and  g  from  Table  III  equation  (38)  yields  the  results  also  shown  in 
Table  III. 


Table 

III,  Values 

and  Results  for  Detector 

Calculations 

Small 

Medium 

Large 

Target 

Target 

Target 

VALUES 

a 

(SR) 

lo"^^ 

-11 

10 

-10 

10 

0 

a 

(radians) 

1.13  X  lO"^ 

3.57  X  lO'^ 

1.13  X  lO"^ 

0 

f 

2.42  X  lO"^ 

3.11  X  lO'® 

6.96  X  lO"^ 

g 

.653 

.842 

1.88 

Results  (a  .-radians 

X  lO"*) 

SI 

M 

=  1 

39.8 

4.17 

0.731 

M 

=  2 

301. 

30.5 

3.45 

M 

=  ^ 

1549 

155. 

16.2 

In  all  cases  the  aperture  was  modeled  as  a  circular  device.  However, 
the  four  sub-aperture  case  (M  =  2)  and  the  12  sub-aperture  case  (M  =  JT?) 
assumed  that  the  circular  aperture  was  subdivided  into  square  sub-apertures. 
The  results  from  Table  III  were  used  directly  for  in  the  single 


aperture  case. 


Now  each  detector  can  measure  only  tilt  (mode  x)  and  tip  (mode  y). 
Thus,  higher  order  modes  must  be  estimated  from  tilt  and  tip  measurements 
made  at  the  various  sub-apertures,  Furthj/r,  if  all  the  sub-aperture 
measurements  have  the  same  noise  behavior,  then  a  least  square  estimate 
is  also  a  minimum  variance  estimate.  This  linear  estimation  procedure 
was  derived  by  seeking  the  sub-aperture  measurements  in  terms  of  the 
first  five  modes  (reference  12).  The  following  measurement  model  was  used. 


m  =  H  a  +  v’ 


(39) 


where 


m  = 


overall  tilt 
overall  tip 
focus 
astig-1 
^astig-2 


tilt(iy 
tilt(2) 

tilt(3)l  and  a 
tilt(4) 
tip(l) 
tip(2) 
tip(3) 
tip(4)  J 

v  =  sub-aperture  noise  vector 
-  m 

m  =  measurements  at  each  sub-aperture 
a  =  optical  modes  into  the  aperture 
The  expressions  for  ^  are 


and 


H.  .  = 
l,j  A 

:l 

•f  ^  dA, 

1  dx 

(40) 

H.  ■  =  -T-  f 

i.J  A.j 

i  Z  ^  =  mode  j 
90®  segments 

f  0^  dA,  4<i<8 

V' 

as  a  function  of 

,  the  H  matrix  is 

X  and  y. 

(41) 

For  a  circle  divided  into 

“i 

0 

2k 

k 

k“ 

1 

0 

-2k 

-k 

k 

1 

0 

-2k 

-k 

-k 

H  = 

1 

0 

2k 

k 

-k 

(42) 

0 

1 

2k 

-k 

k 

0 

1 

2k 

-k 

-k 

0 

1 

-2k 

k 

-k 

0 

1 

-2k 

k 

k 

where  k  =  8/3n 
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Finally,  the  least  square  estimates  are  obtained  from  the  sub-aperture 


measurements  of  tip  and 

tilt  by 

a  =  (h’'  H)'^ 

m 

(43) 

The  covariance  of  the  least  square  estimates  is  (for 

four  identical 

detectors) 

s  = 

(44) 

Thus,  the  variances  for 

the  estimations  are 

(tilt)  =  0^ 

2/4 

(45) 

a2  (tip)  =  0. 

2/4 

(46) 

(focus)  =  o 

i2/32k2) 

(47) 

(astig-1)  = 

0.2/(8k2) 

(48) 

(astig-2)  = 

0.2/(8k2) 

(49) 

where  o.^  =  variance  of 

1 

tilt  and  tip 

for  each  detector.  The  resulting 

sigmas  for  each  mode  and 

the  various 

target  sizes  are 

shown  in  Table 

IV. 

Table  IV. 

Sigma  (x  10  ®)  in 

Radians 

Small 

Medium 

Large 

Target 

Target 

Target 

a(tilt) 

151 

15.2 

1.73 

a(tip) 

151 

15.2 

1.73 

a(focus) 

62.8 

6.34 

.718 

a(astig-l) 

126. 

12.7 

1.44 

a(astig-2) 

126. 

12.7 

1.44 

IX.  RESULTS  AND  CONCLUSIONS: 

Optimal  Control  Results 


The  and  matrices  of  equation  20  and  alpha  of  equation  22 
were  selected  to  produce  a  closed  loop  system  with  a  damped  natural 
frequency  of  about  500  Hz  and  a  damping  ratio  of  about  0.707.  The 
values  used  and  the  resulting  three  controllable  eigenvalues  are  shown 
in  Table  V. 


Table  V.  Values  and  Results  for  Optimal  Control 


Tilt  (&  Tip) 

Focus 

Astig-1 

Astig-2- 

Qjd.i) 

.145x101° 

.2x10° 

.145x101° 

.13x101° 

Ql(5,5) 

121 

95.5 

121 

119 

Q,(6.6) 

. 139x101^ 

.89x101® 

. 139x101’ 

. 135x101° 

. 125x10® 

.6x10® 

. 125x10® 

. 143x10® 

Alpha 

45 

90 

40 

21.1 

Real  Eigenvalue 

-4383 

-7825 

-4367 

-4460 

Complex 

-2215 

-2220 

-2215 

-2222 

Eigenvalues 

±j2218 

±j2220 

±j2218 

±j22l8 

NOTE:  Q^(4,4)  =  Q^(l,l)  =  -Q^(l,4),  -Qj(4,l) 
Eigenvalues  in  per  second 
*For  approximate  3  pole  model 


Kalman  Filter  Results 

The  steady-state  Kalman  filter  was  used  to  estimate  the  states  of 
the  turbulence  model.  As  discussed  in  Section  V  observability  requires 
that  state  (the  integral  of  x^  -  x^)  be  measured.  So,  in  lieu  of 
attempting  to  establish  one  sigma  for  x^,  the  sigma  for  was  varied. 

The  computer  programs  presented  in  reference  14  were  modified  slightly 
to  solve  the  expression  of  equation  (36)  and  to  calculate  a  variance 
reduction  ratio  for  each  set  of  input  data. 

The  results  of  this  research  are  shown  in  the  graphs  of  Figures 
2  through  7.  The  graphs  of  Figures  2  thorugh  4  present  the  results  for 
one  and  four  sub-apertures  and  for  the  three  different  target  sizes 
with  no  time  delay.  The  graphs  of  Figures  5  through  7  present  the  same 
results  for  a  time  delay  of  one  millisecond. 

These  graphs  reveal  that  if  the  variance  of  the  integral  of  the 
error  is  small  enough,  then  four  sub-apertures  provide  a  smaller  normalized 
phase  variance  than  does  a  single  aperture. 

Optical  Detector  Assumptions 

1.  =  4.17  X  lo"** 

2.  A  quarter  circle  can  be  modeled  by  a  square  with  the  same  area. 

(assumed  for  the  four  sub-aperture  case) 
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Control  System  Assumptions 

1.  The  mirror  dynamics  permit  a  damped  natural  frequency  of  500  Hz 
with  a  damping  factor  of  0.707. 

2.  The  modes  are  not  coupled. 

3.  The  time  delay  can  be  modeled  by 

y(t)  =  M  X  (t-T)  +  v^(t-T)  (35) 


15 


Figure  3,  Normalized  Variance  for  Medium  Target,  No  Delay. 


Figure  4.  Normalized  Variance  for  Small  Target,  No  Delay. 


Figure  5.  Normalized  Variance  for  Large  Target,  Delay  =  1  msec. 


MEDIUM  TARGET 

delay  —  1  msec 
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Normalized  Variance  for  Medium  Target.  Delay  =  1  msec. 
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